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Abstract 

The associated production of a Higgs boson with a b quark is a discovery mode for an MSSM 
Higgs boson at large tan /3. We compute the SUSY QCD corrections from gluino and squark loops to 
this process and combine them with the 0{a'^) NLO QCD corrections to obtain reliable predictions 
for the rate. Finally, we compare our results with an effective Lagrangian approximation which 
includes only the low energy corrections from squark and gluino loops to the bb Higgs vertices. 
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I. INTRODUCTION 



The search for the Higgs boson is among the most important challenges facing the current 
generation of colliders. If there is a single Higgs boson with the properties predicted by the 
Standard Model, we expect that it will be discovered at either the Tevatron or the LHC, with 
the optimal channel for discovery depending strongly on the Higgs mass. In the Standard 
Model, the production of a Higgs boson in association with b quarks is never important. 
However, in the minimal supersymmetric model (MSSM), there are five Higgs bosons-two 
neutral Higgs bosons, and H^, a pseudoscalar, A^, and two charged Higgs bosons, H^- 
and the strategy for discovery is quite different from in the Standard Model. In the MSSM, 
the couplings of the Higgs bosons to b quarks can be significantly enhanced and for a large 
range of parameter s]gacej_JIiggs production in association with b quarks is the most likely 
discovery channel 



a fly 

production of a Hi^ 

fiaaflflflfly, 



The production of a Higgs boson in association with a b quark has been extensively 

^^^^^^^^ ^^^^^ ^^^^^ ^ 

studiedlll, 12, 13|, U, 15|, 16|, 17|, 18|, 19|, llO|, luj and the NLO ^CD corrections are well understood. 

In the 4- flavor number scheme. 



both in the 4- and 5- flavor number parton schemes|l|, y, fl. 
the lowest order processes for producing a Higgs boson and a b quark are gg — > bbcp and 
qq — >■ bb(j) ^, fl (The neutral Higgs bosons are generically, = h^, H^, A^). In a 5- flavor 
number scheme, the b quark appears as a parton and potentially large logarithms of the 



form ln(— ) are absorbed into b quark parton distribution functions (PDFs) 
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13|. In the 



5- flavor number scheme, the lowest order process for producing a Higgs boson in association 
with b quarks is bb ^ 6 when no b quarks are tagged in the final state and bq —>■ bch when a 

nnnnn 

single b quark is taggedpj, |3|, 14, |6|, llOj • Contributions from the gg and qq initial states are thus 



subleading in the 5— flavor number scheme. Although the 4- and 5- flavor number schemes 
represent different orderings of perturbation theory, the two schemes have been shown to 
yield equivalent numerical results. The residual renormalization and factorization scale, 
scheme, and PDF uncertainties of the NLO calculations result in a 10 — 20% uncertainty on 
the prediction of the rate[ll, yj. 

When no b quarks are tagged in the flnal state, the dominant production mechanism 
in the 5-flavor number scheme is hb 0, which has been calculated to NNLO ij]. When 
a single b quark is tagged, the rate is lower, but so also is the background and it is this 
channel which we focus on here. Both CDF and DO have derived limits on tan/5 from 
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Tevatron data, based on the search for pp hcj) with cf) t^t~ and cf) — > 66 15|]. In this 
paper, we present the 0{a1) SUSY QCD (SQCD) corrections from gluino-squark loops to 
the h- Higgs production processes. Gluino-squark virtual contributions have been shown to 
be significant for the inclusive bb (p channel 16|] and it is of interest to examine the SQCD 
corrections for the case where a single b quark is tagged along with the Higgs boson. ^ These 



previously computed, ly, and the 



2l[ |. can then be used to obtain reliable 



corrections, along with the NLO QCD corrections 
summation of the logarithmic threshold corrections 
0(q;^) predictions. 

An effective Lagrangian approach has been extensively used in the literature to include 



SQCD effects to the bb(f) vertex 22, l23|, since these corrections are enhanced for large 



values of tan/5. The effective Lagrangian is derived assuming « Msusy, where Msusy 
is a typical squark or gluino mass scale. We investigate the range of validity of the effective 
Lagrangian approach for computing the SQCD corrections to the b(f) production process. 
In addition, we consider the decoupling properties of our results for heavy MSSM mass 
scales. For the decay process, h'^ — > bb, heavy SUSY particles do not decouple unless the 
pseudoscalar mass. Ma, is also large with respect to the electroweak scale 25|, 1261] and we 
observe a similar phenomena in our results. 

In Section 2, we discuss the general MSSM framework of our calculation and the effective 
Lagrangian approach for approximating the SQCD contributions of squarks and gluinos. 
Section 3 contains the lowest order results for the production process, bg b(j), a discussion 
of the renormalization framework for the SQCD NLO contributions, and the SQCD NLO 
results. Numerical results for the Tevatron and the LHC, along with a discussion of the 
decoupling properties of the SQCD contributions are given in Section 4 and analytic results 
for the SQCD corrections are gathered in an appendix. Finally, in Section 5, we present 
some conclusions. Section 5 contains predictions for the b(f) cross sections at the Tevatron 
and the LHC which contain all known Ola"^) contributions and thus are the most reliable 
calculations available. 



^ Correction to Higgs plus jet production in the MSSM from the gg initial state have been examined in 
Refs. 0, 



19 



2C|. 
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II. FRAMEWORK 



A. MSSM Basics 

The MSSM has been extensively studied in the hterature and comprehensive reviews can 



27 



281]. Here we briefly summarize those aspects of the MSSM relevant 



be found in Refs 
for our calculation. 

In the MSSM, there are two SU(2)l Higgs doublets, Hu and Hd, which can be written 



as 



27|,l28|,l29|, 



After spontaneous symmetry breaking, the W and Z bosons obtain masses and fix = 
f ^ + f 2 = (246 GeV)^, while the ratio of the VEVs, tan/5 = ^, is a free parameter of the 
theory. There are five physical Higgs bosons, and , remaining in the theory. 

The scalar potential of the MSSM is described at tree level by two free parameters, which 
are usually taken to be tan/? and the mass of the pseudoscalar Higgs boson. Ma- In terms of 
these parameters, the remaining scalar masses, MH±,Mh and M^, are predicted quantities 
with, 

= Ml + . (2) 
The physical neutral Higgs bosons, and H^, are linear combinations of and h^, 

(3) 

where = cos a and = sin a. The neutral Higgs boson masses at tree level are given by, 

Kh = 1[m1 + MIt 7(Ml + M|)2-4M|Mlcos2 2/5} . (4) 

Eq. m implies a tree level upper bound on the lightest Higgs boson mass, Mh{tree) < Mz- 
Furthermore, at tree-level, 





^1 Ml+ Ml 

tan2a = tan2/j( ^^_^f ). (5) 
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The predictions of Eqs. H] and [5] receive large radiative corrections of 0{GFmf) which raise 
the lower bound on rrih to 130 — 140 GeV.^ We include these corrections using the program 
FeynHiggs, which generates ^ an effective mixing angle, aeff, and radiatively corrected values 
for the Higgs boson masses[30]. Using a^jf in the tree level couplings incorporates the bulk 
of the MSSM corrections to the Higgs masses and mixing angles [3 II]. 
The MSSM Yukawa couplings are given at tree level by, 

Lyuk = -Xbi^iHdbR - XtipiHutn + h.c. , (6) 

where xp]^ = (ti, b^). Eq. [6]generates masses for the t and b quarks, 

XbVi 

= (7) 

We note that the bottom quark only couples to H^, while the top quark only couples to H^. 
In terms of the physical Higgs mass eigenstates. 



(-^bbh° + HbH' - Ubt^^bA'') 



VSM 

^ Hl,?<P^ ' (8) 

where cp = cos/?, S/3 = sin/3, and = tan/5. We see that the Yukawa couplings to the b 
quark are enhanced for large values of tan /3. 

In the MSSM, the scalar partners of the left- and right- handed b quarks, 6^ and bn, are 
not mass eigenstates, but mix according to. 




Lm = -ibhb*j,)M? I r I . (9) 

The b squark mass matrix is, 

' rhl m,(A-Atau;3) ^^^^ 



m{,(yl5 — /itan/5) m 



^2 

R 



2 



The upper bound on the hghtest Higgs mass depends on the stop mass, as weh as other MSSM parameters. 



For a review see Ref. 28 
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where we define, 



;c2 



m 



™2 



MQ + ml + Ml cos 2(3{ll - Qb sin^ 9w) 



m~j^ = Mf^ + ml + Mj cos 2j3Qb sin^ Ow 



(11) 



and Mq^d are the soft SUSY breaking masses, /g = —1/2, and Qb = —1/3. The parameter 
Ab is the trilinear scalar coupling of the soft supersymmetry breaking Lagrangian and fi is 
the Higgsino mass parameter. The b squark mass eigenstates are bi and 62 and define the 
6-squark mixing angle, db 

bi = cos Obbi + sin Obbn 
62 = - sin Obbi + cos 9bbn , 

(12) 



where at tree level, 



tan 29b 



mb{Ab — /itan/3) 



^2 



m 



^2 



m 



R 



Eq. [in] gives the squark mass eigenstates at tree level. 



^ii2 ~ I^^lCos ^b + m^sin 6b] =F sin(26'fe)m6(ylb — /itan/3) 



(13) 



(14) 



The b squarks (unlike the b quark) couple to both Higgs doublets, Hu and and the 
eynman rules for the squark-squark-Higgs couplings can be found in the appendices of Ref. 



271] . Finally, we use the Feynman diagramatic techniques of Ref 3^ and the 



for the squark/gluino interactions which are given in Appendix A of Ref. 33] 



^eynman rules 



B. Effective Lagrangian Approach 



At tree level, there is no ipibuHu coupling in the MSSM (see Eq. [6]), but such a coupling 
arises at one loop through the Feynman diagram shown in Fig. [Hand gives rise to an effective 



interaction 



22 



23j, 



Eq. [15] shifts the b quark mass from its tree level value, 

\bVi _ 



nib 



V2 



;i + Arub) 



(15) 



(16) 



FIG. 1: SUSY QCD contribution from gluinos and squarks to the hhhP effective vertex. Both 
squarks, bi and 62, contribute to the interaction. 



and also implies that the Yukawa couplings of the Higgs bosons to the h quark are shifted 



from the tree level predictions. 



effective Lagrangian approach 22, 



'eff 



le shift of the Yukawa couplings can be included with an 

m 

1 \ / sma 



^^5J\/ V 1 + Arrif, / \ cos/5 
rrib / 1 \ / cos a 



VSM V 1 + 

nib I 1 



cos (5 
— tan (3 



1 + 



tan p tan a J 
Ambtana\-^^o 



1 - 



tan (3 J 
Arrib \y 



vsM V 1 + ^"^f> y V J \ t^'^^ 

= 9bbhbbh' + gbbH^bH' + gbi^lM"" . 



(17) 



The : 
tan P 



jagrangian of Eq. [T7] has been shown to sum all terms of 0(a" tan" (3) for large 
22|.'^ This effective Lagrangian has been used to compute both the inclusive production 



process, bb h'^, and the decay process, — >■ 66, and yields results which are within a few 
percent of the exact one-loop SQCD calculations [l6|, |26|. 

The expression for Am^ is found in the limit ■mh_ << m^ , Mz << mi^,mi^,mg 



The 



contribution to Arrib from sbottom/gluino loops is [221. l23l. |34| 

2as(/ii?) 



Anib -- 

where the function I {a, b, c) is 

1 



37r 



-mgH tan (3 1 {m^-^,m^'^,mg J , 



(18) 



I(a, 6, c) 



(a2-62)(62-c^)(a2-c^)r'^'^°H^^+^'^'^°' 



2 2 1 
c a loe 



(19) 



^ It is also possible to sum the contributions which are proportional to At,, but these terms are less important 
numerically [2^ . 
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6 — r h 



FIG. 2: Feynman diagrams for g{qi) + b{q2) —>■ b{pf,) + (l>i{ph), where 4>i = h^,H^ or A^. 




and asifiR) should be evaluated at a typical squark or gluino mass. Note that Eq. [T8l is 
valid for arbitrary values of tan (3. There is also an electroweak correction to Arrih which we 
neglect, since we are concerned only with the O(a^) contributions in this paper. 

Eq. [18] is a non-decoupling effect in the sense that if the mass scales of the squarks and 
gluino, along with the mixing parameter /i, become large for fixed Ma, Arrib does not vanish, 

Amf, — i> — ^^(^^(/i)— ( tan/? + cot a ) . (20) 
Svr V J 

In the large limit, 

tan/? + cota ^ -^^tan/3cos2/3 + Cj^^J , (21) 



and the decoupling limit of the MSSM is recovered |25[]. 

In the next sections, we investigate the use of the effective Lagrangian of Eq. [17] to 
estimate the SQCD corrections to the process gb —>■ hcj) and compare the results with the 
complete 0{al) one-loop SQCD calculation. 

III. gh 60 AT NLO IN SUSY QCD 

In this section, we summarize the expressions for the SQCD corrections to the gh b(j) 
process. 

A. Lowest Order 

The tree level diagrams for g{qi) + b{q2) —>■ b{pb) + (piph) are shown in Fig. [2] The lowest 
order amplitude can be written as a linear sum over two allowed kinematic structures, with 
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a third contributing at NLO, 



The Born level amplitude is 



where, 



= u{Pb){/li + /l2)l^u{q2) 
Ml = q2uiPb)u{q2) ■ 

"^Born — 9sT gJl/^LO ' 



s t 

and we define the usual Mandelstam variables as: 

s = {qi + q2f = {Pb + Phf ■ 
t = {qi -Pbf = {q2-Phf 

\2 



U 



{Qi - Ph) = {q2 - Pb) ■ 



The spin and color averaged tree level partonic cross section is, 

2 



da 



Born , 



as{,lJ^R) I LO 

9, 



|2 



(22) 

(23) 
(24) 



(25) 



(26) 



dt 96s2 \^bb^ 

The lowest order, 0{as), hadronic cross section is found in the usual manner by integration 
with the b quark and gluon PDFs, 



c^LoiPP ^b(p) = j dxidx2 g{xi, fj.F)b{x2, fJ^p) j + (^1 ^ ^2) 



(27) 



B. NLO QCD Corrections 

The NLO pure QCD 0(«^) corrections to the gb bcf) process consist of one-loop virtual 
corrections containing b quarks and gluons and real gluon emission diagrams, along with 
appropriate counterterms. These corrections have been computed by several groups with 



excellent agreement between the groups [ 
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Schematically, 



(JNLoiPP bcf)) 



QCD 



LO 



OLoipp H) + ^<^qcd{pp b(t>) 



bbfj) ■ 

o'TVLO (gluon only) . 



(28) 



For pseudoscalar production, the three kinematic structures of Eq. include a 75. 
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In Eq. [281 we normalize the Yukawa couplings for both the Born contribution to the NLO 
cross section and for daqco with g^-f^^ (defined by the effective Lagrangian of Eq. [T7j) . 
This algorithm for including the Yukawa couplings of Eq. [T7]in the NLO QCD corrections 
corresponds to the convention of Ref. jl| and so the NLO QCD corrections presented in 
Section 4 (and labelled "NLO (gluon only)") can be directly compared with this reference. 



C. Improved Born Approximation 

One of the major goals of this work is to investigate the accuracy of the Improved Born 
Approximation (IBA). In order to make the comparison as transparent as possible, we define 
an Improved Born Approximation in which the Born amplitude is normalized by the Yukawa 
couplings, g^,l^, of Eq. [13, 

dd'jBA da Born 



dt - dt ' ^^^^ 

The Improved Born Approximation incorporates the effective Lagrangian approximation 
to the SQCD effects on the 660 Yukawa couplings at low energy, but does not include 
the full SQCD calculation presented in the next section. In particular, the "Improved Born 
Approximation" does not include contributions from box diagrams including internal squarks 
and gluinos or the full momentum dependence of the SQCD contributions. 



D. One Loop SQCD Corrections 



The amplitudes for the SQCD one-loop corrections from sbottom/gluino loops can be 
written as: 



(30) 



where the matrix elements Ms,Mt and Mi are given in Eq. [22] and expressions for the 
individual contributions, Xj are found in Appendix B. The integrals of Appendix A are 
evaluated in n = 4 — 2e dimensions. The diagrams consist of self-energy, vertex, and box 
contributions, along with the appropriate counterterms to cancel the ultraviolet divergences. 
We consistently neglect contributions which are suppressed by terms of 0(m^/s, m^/M|), 
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etc. The spin and color averaged partonic cross section to NLO in SQCD is, 



xf'M; + X\"'Ml ] + 



(1), 



da, 



CT 



dt ' 



(31) 



where the contribution from the counterterms is discussed in the next section. The SQCD 
NLO contribution to the hadronic cross section is 



O-NLoiPP H)SQCD = / dXidX2 



g{xi,^)b{x2,fi) 



hi 



(32) 



The resuh of Eq. [32] is labelled "NLO (gluino/squark only)" in the figures of Section IV. 
Finally, we note that 

'tmj, + + us' 



\Alo-m:\ 

I^lo-m;| 

\Alo-M*\ 



-4 
-4 
-2 



t 

sw?, — us -\- umi 



u + us 
t 



(33) 



E. Counterterms 



In this section, we discuss the squark and gluino contributions to the counterterms which 
are necessary for the SQCD one-loop calculation. The renormalization of the pure QCD 
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contribution has been presented previously [l|, I^, 0, [s, 7, 

Self-energy and vertex corrections to the tree-level process gb bcf) process give rise 
to ultraviolet (UV) divergences. These singularities are cancelled by a set of counterterms 
fixed by well-defined renormalization conditions. The renormalization of the propagators and 
interaction vertices of the theory reduces to introducing counterterms for the external field 
wavefunctions of the bottom quarks {5Zv) and gluons (^Zs), and for the strong coupling 
constant {5Zg). The counterterm for the bottom quark Yukawa coupling, g^-^ ~ Tn^/v, 
coincides with the counterterm for the bottom quark mass, since the vacuum expectation 
value V is not renormalized at the one-loop level of SUSY QCD. (We follow the approach of 



Refs. 



33 



35] for the SQCD renormalization). 
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b 



9^ 



FIG. 3: SUSY QCD contribution to the b quark self energy. Both squarks, bi and 62, contribute 
to the interaction. 

We define the b quark self energy as 



which yields the renormahzed propagator, S^g„(p) = Y^'ip) + 



6mt 
nib 



(34) 



(35) 



We evaluate the SQCD contribution to the b mass using the on-shell renormalization con- 
dition, 



mb) = 



lim. 



i)-mb 



0. 



(36) 



This scheme decouples the gluino and 6-squark from the running of the b quark Yukawa 
coupling. The self energy contribution from gluino-squark loops is shown in Fig. [3] and 
yields the on-shell result. 



rrib 



a, 



Stt t 



+2 sin 9bCos9b 



nib 



Bo{nil,nig,nij^J - Boinil,ni-g,nij^^] 



(37) 



where the integrals Bo and Bi are defined in Appendix A. 



The 6-quark self-energy is renormahzed by \E' — > a/ Zy"^ = ^1 + 5Zy'^. Neglecting 
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contributions suppressed by powers of mj,, in the on-shell scheme we find, 
6Zt = -S 



V ~ \p^=ml 

2 ™ r ^ , D 



37r \ ^ 



(38) 



The external gluon has the wavefunction renormahzation, \fZlg^ = y/1 + SZ^g^ 

and the renormahzation of the strong coupling is found from g^ —>■ Zgg^. We have fur- 



ther, 5Zg = —5Z3/2} We renormalize gs using the MS scheme modified to decouple 
heavy SUSY particles, i.e. the heavy squark and gluino contributions are evaluated at zero 



momentum 36] 



SZ. 



\pole 



dp 



2 



where Ej^ represents the transverse portion of the gluon self energy, the sum in Eq. [39] is 
over all squarks (not just the b squark), and fiR is an appropriate renormahzation scale. 

The Yukawa couplings are defined by EqlTTl which includes the summation of large tan f3 
effects. Using Eq. [T7| for the Yukawa couplings of the Born contribution in Eq. [31] includes 
some one-loop effects in the first term which must be subtracted in order not to double 
count. This generates additional counterterms, 

drhb'' ^ A 1 



Ami, ( 1 H (h production) 

nib V ^ P J 

Sfhb^ A A tana\ , , . , 

= Anib 1 [rl production) 

nib V P / 

The counterterms make a contribution to the total partonic SQCD cross section, 

^ = ^ f 26Z^y + 6Z, + 26Z, + 2^ + 2^) . (41) 
at at \ nib nib J 

We note that the 1/e contributions cancel in Eq. [HI although there are remaining finite 

pieces. 



5 



The validity of this relation for the squark/gluino contributions to the SQCD renormahzation is demon- 
strated in Ref. 



33[ and references therein. 
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F. Complete NLO Result 



The complete Ola"^) NLO prediction includes both the pure QCD NLO corrections and 
the squark/gluino NLO SQCD contributions, along with the Born contribution® , 

(^NLoiPP H)qCD+SQD = (^NLoiPP ^ H)sQCD + +^(yQCD ■ (42) 

The curves labelled "Complete NLO" in Section IV correspond to Eq. |32] and represent the 
most accurate predictions available for the pp —* hcf) MSSM cross sections. 



IV. RESULTS AND DISCUSSION 



Our NLO numerical results are obtained with CTEQ6M parton distributions [37|], as^Hn) 
evaluated with the 2— loop evolution and a^^^{Mz) = 0.118, and use the Yukawa couplings 
of Eq. [T71 The lowest order cross sections use CTEQ6L PDFs, asifiR) evaluated at 1-loop, 
and have the lowest order Yukawa couplings of Eq. [HI We require that the outgoing b quark 
has pt > 20 GeV and pseudorapidity \ f] \< 2.0 for the Tevatron and | |< 2.5 for the LHC 
In the NLO QCD real gluon emission contribution (obtained from Ref. (6[), the final state 
gluons and b quarks are merged into a single jet if their pseudorapidity /azimuthal separation, 



AR = a/ (At])^ + (A0)^, is less than 0.4. The renormalization/factorization scales, ^r,^Fi 
are taken to be Finally, the results labelled "Improved Born Approximation" (Eq. 

[291) . use NLO PDFS and asifiRj evaluated with the 2— loop evolution and a^^^{Mz) = 



0.118. 



The b quark mass appearing in the Yukawa couplings of Eqs. [S]and[T7|is taken to be the 



running MS b quark mass and is evaluated at two loops for the NLO predictions and the 
IBA predictions and at 1— loop for the LO predictions, 



m 



pole 



m 



pole 







1 pole\ 






co/6o - 


1 pole\ 





1 + 



Co 



ci - bi 



TT 



(43) 



Note the the Born contribution is (arbitrarily for this purpose) assigned to ctnlo^P 



)SQCD 
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TABLE I: b squark masses and mixing angles from 30|]. All soft-SUSY breaking masses are taken 



A), = At and /i = 


= M2 = 


200 GeV. 






tan P 


= 10 




MsusY (TeY) 


h m-,^ (TeV) 




1 


-0.64 


1.0 


1.0 


2 


0.00 


2.0 


2.0 


3 


0.64 


3.0 


3.0 


4 


0.71 


4.0 


4.0 


5 


0.74 


5.0 


5.0 



TABLE IL b squark masses and mixing angles from 30]. All soft-SUSY breaking masses are taken 



equal to MsusY = w-s = -A-b = -A-t and /U = M2 = 200 GeV. 



tan /3 = 40 



M5c/5y(TeV) h m-,^ (TeV) m^^ (TeV) 



1 


-0.76 


1.0 


.98 


2 


-0.76 


2.0 


2.0 


3 


-0.75 


3.0 


3.0 


4 


-0.75 


4.0 


4.0 


5 


-0.73 


5.0 


5.0 



where, 



co = i, ci = ^(mN,-10Nf, 



with Nc = 3 and Nfi = 5, the number of light flavors. We take the b quark pole mass, 
^poie ^ 4_g2 GeV. 

The MSSM parameters are found using FeynHiggs to generate an effective Higgs mixing 
angle, a^ff, and radiatively corrected Higgs masses. The MSSM parameters are listed in 
the figure captions. The gluino mass is always identified with Msusy, while the b squark 
masses and the b squark mixing angle, 6b, are given for representative parameter values in 
Tables [T] and [m 
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FIG. 4: Scale dependence of the NLO results for pp hhP (LHS) and pp bH^ (RHS) for 
tan/3 = 40. We set the renormalization/factorization scales equal, fin = fip. 

The NLO results for the pp{pp) —>■ hcj) processes depend on both the renormalization 
and factorization scales. For simplicity, we take these scales equal, /ij? = ^f- Fig. H] 
shows the dependence on scale (/ii?//io) of the NLO results at both the Tevatron and the 
LHC and demonstrates the expected improvements between the lowest order results and the 
NLO results. The curves labelled "Complete NLO" include the full set of QCD and SQCD 
contributions (Eq. I42l) and are significantly less dependent on the choice of scale than the 
lowest order result (Eq. [27]). The dashed curve ("NLO (gluon only)", Eq. [28]) only includes 
the SQCD contributions via the effective Yukawa couplings of Eq. [T7] and has an almost 
identical scale dependence as the full result of Eq. [121 For the remainder of our plots, we 
choose the renormalization/factorization scale to be 

Results for the lightest MSSM scalar, h^, are shown in Figs. [5] and [6] for the Tevatron.^ 
Fig. [5] shows the dependence of the complete set of SQCD contributions (Eq. [32]) and the 

^ A preliminary study of the SQCD contributions to 6/i° production in the MSSM appears in Ref. (ssj . 
This work used a renormahzation framework in which the heavy squarks and gluino do not decouple and 
hence their results are not directly comparable to ours. We find, however, qualitative agreement with 
their results. 
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pure NLO QCD result (Eq. [28]) on the lightest Higgs mass for tan/? = 10 and 40. The curve 
marked "gluon (only)" (Eq. [281) includes the NLO QCD corrections and uses the effective 
Yukawa couplings of Eq. [TTl while the SQCD curve (Eq. [521) incorporates the full set of 
SQCD corrections (box diagrams, etc). We see that the effects of squark and gluino loops 
which are not absorbed into the the effective Yukawa couplings are typically less than 1 — 2% 
for both tan (3 = 10 and tan /? = 40 and can safely be neglected. We also note that the sizes 
of the NLO corrections (both pure QCD and SQCD) do not have strong dependences on 
the mass of the produced Higgs boson. 

Fig. [6] compares the SQCD result (Eq. [321) with that obtained in the Improved Born 
Approximation (IBA), Eq. [221 For both values of tan/3 plotted, the Improved Born Ap- 
proximation is an excellent approximation. This is an important result of our calculation 
because it demonstrates that the effective Lagrangian approach is extremely accurate for 60 
production. It is apparent that the squark and gluino contributions which are not included 
in the IBA approximation give a negligible contribution to the rate. Figure [6] also exhibits 



the slow approach to the decoupling limit already noted in Ref. [251]. The decoupling limit, 
where the SQCD corrections do not contribute to the rate, corresponds to o'nlo/'^lo — ^ = 
in this plot. The approach to the decoupling limit is significantly slower for large tan/5 than 
for small tan f3, as is apparent in Eq. [201 Even for Msusy as large as 5 TeV, the SQCD 
effects present in the couplings of the effective Lagrangian of Eq. [17] are still (9(5%) for 
tan (3 = 40. 

Results for the heavier neutral MSSM scalar, H^, are shown in Figs. [7] and [HI for the 
LHC. Fig. [7] shows the dependence of the SQCD corrections on the heavy Higgs mass for 
tan/3 = 10 and 40. The NLO pure QCD results (Eq. [28]) are quite sensitive to M^, while 
the SQCD corrections (Eq. [32]) are relatively independent of the Higgs mass and quite small. 
The complete NLO rate including all QCD and SQCD effects (Eq. [^2]) differs by less than 
1 % from the "NLO (gluon only)") calculation where the SQCD effects are included only 
via the Yukawa couplings of the effective Lagrangian. Fig. [8] compares the full SQCD result 
with that obtained in the Improved Born Approximation (IBA) and as is the case for the 
lighter Higgs boson, the Improved Born Approximation closely approximates the full SQCD 
calculation. The slow approach to the decoupling limit for the is similar to that seen for 
the h^. 

Results for the MSSM pseudoscalar, A'^, are shown in Figs. [91 and [TOl for the LHC. Fig. 
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FIG. 5: Effects of NLO QCD and SQCD contributions on the rate for pp bh^ at the LHC. The 
outgoing b quark satisfies pT > 20 GeV and | r/ |< 2. The NLO results use NLO PDFs, 2— loop 
evolution of as{fJ-R) and fnh{fiji), and the Yukawa couplings of Eq. [TTl The plots are normalized 
to the lowest order cross section of Eq. [26l which is computed with lowest order PDFs, 1-loop 
evolution of as{nR) and rrTh{fiji), and the lowest order Yukawa couplings, 5^^- 

|9] shows the dependence of the SQCD corrections on the pseudoscalar mass for tan /5 = 10 
and 40. The SQCD effects (Eq. [321) are positive and are of a similar magnitude as those 
found for production and again are relatively independent of the Higgs mass, Ma- Fig. 
[To] compares the full SQCD result with that obtained in the Improved Born Approximation. 
In this case, the Improved Born Approximation gives a prediction which differs from the full 
SQCD rate by about —3% for both tan/5 = 10 and 40 and we again see the slow approach 
to the decoupling limit. This plot fixes = 500 GeV and varies Msusy- From Eq. [211 it 
is clear that for decoupling to occur we need to also take — > c>o. For fixed Ma, even for 
Msusy = 5 TeV, the effect of squark and gluino loops are still 0(5 — 10)%. 
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FIG. 6: Comparison of the full SQCD calculation (Eq. [32]) with the Improved Born Approximation 
(Eq. [29l) on the rate for pp hhP at the Tevatron. The outgoing b quark satisfies pt > 20 GeV 
and I ?7 |< 2. 

V. CONCLUSIONS 

We have computed the 0{al) SUSY QCD corrections from squark and gluino loops for 
the associated production of an MSSM Higgs boson and a b quark. When the Yukawa 
couphngs are normaUzed with the effective Lagrangian of Eq. [TTl the remaining SQCD cor- 
rections from squark and gluino loops are typically of order a few percent. We have therefore 
explicitly demonstrated that the effective Lagrangian approach to calculating the effects of 
squark and gluino loops is extremely accurate and can be reliably used to approximate the 
b(f) cross sections. In Figs. [H] and [121 we summarize our results by showing the complete 
NLO predictions for b Higgs associated production at the Tevatron and the LHC^ These 
figures include all NLO O^a"^) contributions and represent the most complete calculation 
available. 

® The bh cross section decreases rapidly as the maximum allowed value of is approached from below due 
to the suppression of the bbh coupling. Similarly, the bH cross section decreases as M^ Max is approached 
from above due to the suppression in this region of the bbH coupling. 
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FIG. 7: Effects of NLO QCD and SQCD contributions on the rate for pp bH'^ at the LHC. The 
outgoing b quark satisfies pT > 20 GeV and | ?7 |< 2.5. The NLO results use NLO PDFs, 2— loop 
evolution of as{fJ-R) and frTh{fifi), and the Yukawa couplings of Eq. [TTl The plots are normalized 
to the lowest order cross section of Eq. [26l which is computed with lowest order PDFs, 1-loop 
evolution of as{nR) and rrTh{fiji), and the lowest order Yukawa couplings, 5^^- 

Acknowledgments 

The authors would like to thank C. Kao and Y. Wang for useful discussions. This work 
is supported by the U.S. Department of Energy under grant DE-AC02-98CH10886. S.D. 
thanks the SLAC theory group for their hospitality, where this work was completed. 



20 



Vs= 14TeV 



tanp = 10 




3000 
MgusY [GeV] 



FIG. 8: Comparison of the full SQCD calculation (Eq. [32]) with the Improved Born Approximation 
(Eq. [29]) for the rate for pp hH^ at the LHC. The outgoing b quark satisfies pT > 20 GeV and 
I r] |< 2.5. 

APPENDIX A: SCALAR INTEGRALS AND TENSOR COEFFICIENTS 



The scalar integrals are defined as: 



Co{pi,pi;Mo,M^,M2) 



Do{pi,pipi;Mo,M,,M2,M^ 



d^k 1 

(P^k 1 

(27r)" A^oA^i ' 
1 

1 

(27r)" AroA^iA^2iV3 



(Al) 



where: 



Ml 



{k + p{)^-Ml 

{k + pi+p2f -Ml 

{k+p^+P2+p^f - Ml 



(A2) 
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FIG. 9: Effects of NLO QCD and SQCD contributions on the rate for pp bA° at the LHC. The 
outgoing b quark satisfies pT > 20 GeV and | ?7 |< 2.5. The NLO results use NLO PDFs, 2— loop 
evolution of as{fJ-R) and frTh{fifi), and the Yukawa couplings of Eq. [TTl The plots are normalized 
to the lowest order cross section of Eq. [26l which is computed with lowest order PDFs, 1-loop 
evolution of as{nR) and rrTh{fiji), and the lowest order Yukawa couplings, 5^?. 



The tensor integrals encountered are expanded in terms of the external momenta Pi and 
the metric tensor g^'^. For the two-point function we write: 

;B'^{pI;Mo,M,: 



IQ^2 v^i' J (27r)«iVoAri 

= ^p'^B,{pIMo,M,), (A3) 

while for the three-point functions we have both rank-one and rank-two tensor integrals 
which we expand as: 

C^\pIpI-Mo,M,,M2) = p'^p\C2i+p^,p',C22 

+ {p'lP'2+PlP'2)C23+9'''C2,, (A4) 
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FIG. 10: Comparison of the full SQCD calculation (Eq. [32l) with the Improved Born Approximation 
(Eq. [29]) for hJ^ at the LHC. The outgoing h quark satisfies -pT > 20 GeV and | r? |< 2.5. 
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FIG. 11: Complete NLO 0{al) results for pp bh^{H^) at the Tevatron. 



where: 



-C7^(Cn(p?,P2;Mo,Mi,M2) 



(A5) 



Finally, for the box diagrams, we encounter only rank-one tensor integrals which are 
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FIG. 12: Complete NLO ©(a^) results for pp hh^{H^) at the LHC. 



written in terms of the Passarino-Veltmann coefficients as: 



D^^{pi,pi;Mo,MuM2) 



(27r)« NoN,N2N3 
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APPENDIX B: ONE LOOP SQCD COEFFICIENTS 



The cofficients for each diagram are given below in the notation of Eq. 
The self-energy contributions are shown in Fig. 




b '\ 
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FIG. 13: Self- energy diagrams, Si and 82- 



Diagram 5*1 : 



24 



Diagram 5*2: 



4f 



4? 

4 
~3 



Bi{t;mg,mf^) + Bi{t\mg,m^;J 



fBll 



X 



is) 
S2 



X 



it) 

52 



X 



^9bbij) 
(1) 

52 







(B2) 



The virtual diagrams are sliown in Figs. [TU [T5] and [ini An effective bbg vertex can be 
extracted from Figs. [TU and [I5] and in tlie limit where the b quarks of the effective vertex 
are on-shell, our results agree with Ref. 33(]. An effective vertex can be found from 
Figs. [16] and agrees with the results of Refs. 22, 23, [25] when the b quarks are on shell. 
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Diagram Vi : 



FIG. 14: Vertex diagrams, Vi and V2 
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(0,0;m5,mj,-,m^-). 



-C24 + (^12 + C23 
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(B3) 



Diagram V2 : 
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= -3-2^— (s.0;n>5,mj,,mj,) 
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Cll — Ci2 + (^21 ~ (^23 



(s, 0; m 
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FIG. 15: Vertex diagrams, V3 and V4. 
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Diagram V3 
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Diagram V4 : 



FIG. 16: Vertex diagrams, V5 and Vg 
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Diagram V5 : 







X 



^('^)<!2cOS^,sill^;. 



+(cos^ ^6 - sin^ 9b)Cii2 
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61' "^61^ 
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(B7) 
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where the full arguments of the scalar integral and tensor coefficients are 
{pj,pl;mlmlml) = (0, t; m|, m?, m|). 
Diagram Vg : 



X 



+ (cos^ 6b - sin^ Ob)Cii2 



Co (m| , m| ) + Co (m| , m| ) 



(1) 



0, 



(B8) 



where the full arguments of the scalar integral and tensor coefficients are 
{pl,pl;ml,ml,ml) = (0, s; , m?, ). 



1. The Box Diagrams 




FIG. 17: Box diagrams Bi, B2 and B-^. 



The box diagrams are shown in Fig. [T71 The results of Eqs. IB9t IBIOI and IBllI sum over 
the contributions of bi and 62- 
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Diagram Bii 
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(B9) 



where the full arguments of the scalar integral and tensor coefficients are 
(pl,pl,ph mi, 1712, 1713,1714) = (0,0,0;mft-.,mg,m|,m|). 
Diagram D2: 
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{Du - D,s) {ml ,m}) + {D,2 - ^^13) K , mj- ) 



(BIO) 



where the arguments of the scalar integral and tensor coefficients are 
(pl,pl,phml,ml,ml,ml) = {0,0,0;mg.,mg.,mg,mr). 
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Diagram B3: 
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where the arguments of the scalar integral and tensor coefficients are 
{Pl,pl,phmi,m2,m3,m4) = {0,0,0;mg.,mg.,mg,mr). 
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